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U: 

• . Abstract 

O ■ Lie group theory states that knowledge of a m-parameters solvable group of symmetries of a system 

of ordinary differential equations allows to reduce by m the number of equation. 

, , We apply this principle by finding dilatations and translations that are Lie point symmetries of 

►^ ■ considered ordinary differential system. By rewriting original problem in an invariant coordinates set 

f~^ ' for these symmetries, one can reduce the involved number of parameters. This process is classically call 

\^ , nondimensionalisation. 

f^ ' We present an algorithm based on this standpoint and show that its arithmetic complexity is poly- 

^^ ' nomial in input's size. 

o: 

O ; 1 Introduction 

^ ' This paper is devoted to the process of nondimensionalisation that it described in ^U] as follow: "Before 

^ , analysing [a] model it is essential, or rather obligatory to express it in nondimensional terms. This has 

several advantages. For example, the units used in the analysis are then unimportant and the adjectives 
small and large have a definite relative meaning. It also always reduces the number of relevant parameters 



X 



C^ ' to dimensionless groupings which determine the dynamics". 

Example 1 — In order to illustrate this statement, let us consider the following Verhulst's logistic growth 
model with linear predation (see § 1.1 in \W\). 

dx/dt — x{a — bx) — ex, a = b = c = 0, (1) 

for which all forthcoming computation could be easily performed by hand. 

Step 1. One can remark that the following one parameter group of translation symmetries: 

Ta : a ^ a + A, c ^ c + A (2) 

and the following two parameters group of scale symmetries: 



a - 


^ va, 


b - 


-> vb/ fj. 


c - 


-> vc, 



leave invariant solutions of system (|T]l. 

Step 2. Assuming that a ^ c^ one can determine some rational invariants of these groups like: 

t = (a - c) i, y = X. (4) 

These invariants are the dimensionless groupings evoked in above citation. 

Step 3. Using above rational invariants, system (^ can be rewritten in these new coordinates as follow: 

dj:/dt = y(l-y). (5) 

In fact, one can deduce from ^ and Q a differential system that could be simplified in order to obtain ^. 
Thus, groups |(5J and (PJ allow to rewrite our original problem with a reduced number of parameters. 

We have considered a variety of such dynamical systems that arise from biological, physical, etc. models. 
As often observed, those systems can be rewritten in a lower number of variables or parameters than initially 
appears. 

There is no difficulty to perform manipulations done above by hand on simple systems. But while 
system's complexity increases, these manipulations require more avoidable work as illustrated by the following 
example. 

Example 2 — Let us consider the following prey- predator model taken from page 88 in |ll)| : 

n. 




(6) 

= 0. 

Step 1. One can determine that the following 3 parameters group of scale symmetries: 

t — > i/A, r -^ Xr, h -^ vh, 
S(\n,v) ■■ n ^ nn, s -^ As, ki -^ ^iki, (7) 

p — > ^J'■p/v, e — > /xe, k2 -^ \vk2 

leaves invariant solutions of system ©. Computations show that this system does not have any translation 

symmetries. 

Step 2 & 3. If group's parameters are specialized as follow: 

^ = ? ^=i' "=^' (^) 

resulting transformation S(i/r^i/ki,rlk2) induces a change of coordinates that send coordinates r, fci and fc2 
to 1 and other coordinates to rational invariants: 

n koV s e , rh ,„, 

fci fcir r ki fc2 

After applying transformation iS(i/r-,i/fci,r/fe2) on system ^, we obtain a new system expressed in above new 
coordinates: 

s = e = f) = 0. 

and as in example ^ we reduce by 3 the number or parameters. Thus, any further manipulations of 
system (^ (for example, its phase plane analysis or study of its bifurcation analysis) is simplified in the new 
set of coordinates ©. 



Same process could be achieved in another set of invariants coordinates by choosing to send s, e and h 
to 1 using transformation S(^i/s,i/e,i/h)- 

Remark 1 — Example |2] shows how scale symmetries allow computation of a set of invariant coordinates 
and of resulting new system by an evaluation without any further algebraic manipulation (see Section [4.2|l . 

1.1 Dimensional Analysis and some Lie point Symmetries 

There is just seven primary units in the International Metric System and all physically meaningful equations 
could be written in these units. Such physical primary dimension (say time) can be expressed in different 
units (second, hour, etc.) and thus, these units can be scaled independently of each other (a hour is 3600 
seconds, etc). This possibility induces also dilatations on secondary dimension like speed [m/ s) or corporal 
mass index (kg/m?), etc. Thus, there is likely some scale transformations group acting on functional relation 
among these kind of quantities. 

Dimensional analysis is based on this remark. Bridgman explains in 3 : "The principal use of dimensional 
analysis is to deduce from a study of the dimensions of the variables in any physical system certain limitations 
on the form of any possible relationship between those variables". Thus, dimensions analysis addresses very 
general problems and we consider in this paper just one of its consequences that is nondimensionalisation. 
This process is done because reducing the number of parameters, or possibly also the number of state 
variables, is an advantage for studying qualitative features of the model. 

Previous related works. For large systems, nondimensionalisation could become a difficult process that 
motivate several implementations (see |5] and the references therein for more details). Up to our knowledge, 
there is no complexity result concerning these works that are related to 11 theorem and to rules of thumbs 
based on the knowledge of units in which is expressed the problem. 

As we notice at beginning of this section, there is often some scale symmetries group of ordinary differential 
system describing biological, physical, etc. phenomena. This fact could be considered in the framework of 
Lie symmetries group theory (see § 3.4 in |11| for such presentation of 11 theorem). This standpoint allows 
to presents computation involved in nondimensionalisation in a very simple and efficient way that we did 
not found in literature. 

The present study was motivated by this fact and by authors' inability to apply classical dimensional 
analysis' rules to systems composed of more then 7 equations and a dozen parameters. 

Remark 2 — In this paper, we restrict ourself to translation and scale symmetries that occurs frequently 
in application. Same type of result could be obtain for other kind of Lie point symmetries (rotation, inversion, 
etc. See Example |31 in Section Isj. 

1.2 Main contribution 

It is the aim of this paper to provide an algorithm to make the reduction evoked in previous section. In 
the example above and 90% that we tried, the system (|ll|l is actually symmetric under a group of scalings 
and/or translations. New variables that we introduced in above examples and many others actually form a 
generating set of rational invariants for some group action. 

Computing symmetry of a differential system has been an intensive field of application of computer alge- 
bra, especially to mathematical physics lllj and computing a generating set of (differential) invariants 
for general group action has received relatively recently firm foundations 7, 6, 11 . 

The viewpoint of this paper is not to be general nor theoretical but practical for a large class of problems. 
We shall apply known general theory to special cases. We obtain efficient algorithms for reducing the number 
of parameters in biological, chemical, etc. models basing ourselves on the observation of a general scenario: 
the invariance of the model under a group of scaling and translation. We provide efficient computer algebra 
algorithms for computing the scaling and translation symmetry of a differential system (|ll(l , compute their 



invariants and rewrite the system in terms of those. We thus obtain the reduced system. These result are 
summarized in the following statement: 

Theorem 1 — Let E be a differential system bearing on n state variables and depending on £ parameters 
that is coded by a straight-line program of size L (see Section |21l. There exists a probabilistic algorithm that 
determines if a Lie point symmetries group of E composed of dilatation and translation exists; in that case, 
a rational set of invariant coordinates is computed and E is rewrite in this set with a reduced number of 
parameters. 

The arithmetic complexity of this algorithm is bounded by 



0[{n + t+l){L+{n + l+ l)(2n + ^ + 1))) . 



Outline of the paper. In the next section, we recall some basic definitions of differential algebra and 
we present in this framework the relationship between dilatation/translation transformations and induced 
derivations that are theirs infinitesimal generators. Then, we recall that such infinitesimal generators are 
defined by a partial differential equations system i.e. some infinitesimal conditions which are presented and 
used in the sequel. 

In the last part of this paper, we show that infinitesimal conditions allow by a Gaussian elimination 
performed on a constant field to determine Lie groups of scale/translation symmetries. In fact, for these 
transformations, infinitesimal conditions split into a linear system of algebraic equations. We point out that 
computation of sets of rational invariant coordinates and rewriting of original system in these sets is reduced 
to linear algebra and we estimate the related arithmetic complexity. 

2 Mathematical framework 

Hereafter, we consider an algebraic ordinary differential system bearing on n state variables X := (xi, . . . , Xn) 
and depending on i parameters :— {0i, . . . ,6i): 



i = i, e = o, 

X ^F{t,X,e). 



(11) 



The letter X stands for first order derivatives of state variables (ii, . . . , in) w.r.t. time t; we use the standard 
notation X'*) = {xi^^^\ . . . , a;„^*^) for higher derivatives of order i. We assume that F := {fi, . . . , fn) consist 
of rational functions over a subfield K (Q for example) of C i.e. F is a finite subset of K(i, X, 8). 

Complexity model. We evaluate the complexity of our algorithms within the model of straight-line 
program (see § 4 in ^). For instance the expression e :— {x + 1)^ is represented by the following kind of 
instructions sequence: 

ei:=a; + l, 62 := ei^, 63 := 62^, 6:^6361. (12) 

The complexity is measured in terms of the following date of the input: n + i and the number L of arithmetic 
operations needed to compute the numerators and denominators of F. 

2.1 Differential Algebraic Setting 

We use differential algebra, founded by J.F. Ritt, in order to introduce forthcoming definitions (see [T5| 
for a complete description). 

The differential algebra K{t,X,Q} is the IK-algebra of multivariate polynomials defined by the infinite 
set of indeterminates {i, 0,X(^^| Vj G N} and equipped with time derivation d/dt denoted by C. Thus, for 
any y in K{t,X,Q}, relations Cy'^^^ — y(*+i) hold. 

System (|ll|l defines a prime differential ideal / of the algebra K{t, X, 9}. This ideal encodes all relations 
between coefficients of power series solutions of E. In the sequel, we are going to focus our attention on the 
quotient differential fraction field K{t, X, 6}// — denoted by K. — associated to /. 



Derivations vector space. Let us recall that derivations acting on K(t, X, <d) form a vector space 
over K(i,X, 0) denoted by Der(K(t, X, 9)/IK) and equipped with a canonical base given by elementary 
derivations: 

This vector space equipped with canonical Lie bracket forms a solvable Lie algebra (see ^1 E] and references 
therein for some algorithmic tools used in study of Lie algebra in this context). In the sequel, we are going 
to consider other such algebras (see remark IT!^ . 

Canonical field isomorphism. As our input system defines explicitly a vector field, any high order deriva- 
tives could be rewritten to 0th order ones using relations Hll|l. Thus, our dilTerential field /C is isomorphic 
to the difi'erential field K(t, X, 8) equipped with the following formal Lie derivation in Der(K(i, X, 0)/K): 

■i— 1 

We are going to use directly the canonical isomorphism between differential fields (A^, C) and (IK(i, X, <d),T>). 
All forthcoming developments are based on derivation V. Thus, hereafter, we denote by K, the field K(i, X, O), 
the set {Vfi, . . . , Vfn) by VF and the composition T) o ■ • -oT) by VK 

j times 

Formal power series. Let us denote by S(t, X, 9) formal power series with coefficients in IC that are 
solutions of S = F{t, S, 8) with initial condition S(0, X, 8) := X. 

These power series could be define using derivation V by the following formal relations: 

5(i,X,8)=^P^X^. (15) 

Remark 3 — Higher order derivatives could be considered via differential field K(t, X, 8) that is the 
field K(t, 8, X'-'^\i G N) and the formal derivation: 

v^=v+ y: e^'/'9^' (16) 

jeN\{o,i} J=i * 

as VooX^ C /C, that allows to compute higher order derivatives of X. 

In the sequel, we are going to exploit the fact that the differential field {IC, V) encodes all formal infor- 
mations associated to formal power series solution of S. 

3 Infinitesimal Generators of Scale Symmetries acting on K 

We are looking for a 77i-parameters group <Jix^^...^\^) acting on /C that is a symmetries group leaving the 
solution of E invariant. First, we are going to recall how a derivation in Der(/C/IK) could be associated to 
a scale (resp. translation) transformation. Then, we explicit that such a derivation defines a scale (resp. 
translation) symmetries of E. 



3.1 Scale/translation transformation groups and associated derivations vector 
spaces. 

In order to explicit classical basis of forthcoming computations, let us consider the following one-parameter 
group of scale transformations: 

f — > A"* f 

xi — > A^i xi, 
^a: . : (17) 

the group parameter is denoted by A and taken in a field F; exponents a^ are in a constant field (Q for 
example) and they define considered group. 

Remark 4 — For translation, we consider analogous expressions y -^ y + ayX with y in (t,X,Q). In 
Example ^ we already considered the one-parameter translation given by the exponents set defined by 
relations aa = ac = ^ and at = ax = at = that leads to transformations 

a\ :t ^' t, X -^ X, a ^ a + A, b ^ b, c ^ c + A. 

Actions of these groups could be defined on power series S solution of E (i.e. they remain invariant under 
group action 'E.{a\{t, X, 6)) = E{t, X, 9)) and on K(t, X, 8). We consider this last case in the sequel. 

Infinitesimal generators vector space. When a transformation group is connected — as it is always 
the case for transformations considered in this paper — each symmetries of group (|17|l is associated to an 
infinitesimal generator (for geometric description, see § 2 in jll|'l. It is a derivation on /C in the vector 
space S generated in Der(/C/K) by the derivation: 

^5= E ^yy^' (18) 

with the ay are constant exponents taken from IjlTI) . In fact, given any element Sp in S equal to pS with p 
an element of a constant field, one can construct a scale transformation using exponential map as follow: 

ax: IC ^ /C(A) (19) 

if exponential A := exp(p) is defined. At opposite, one can determine a derivation associated to the applica- 
tion (7\ but we are not going to use this fact in the sequel. 

Remarlc 5 — In the case of a m-parameters group, the vector space S is generated by m deriva- 
tions Si, ... , Sm of the same type then H18|l . More precisely S is a Lie algebra s.t. [5^, tSj] = 0. 

Remarlc 6 — Same considerations are true for translation but in this case infinitesimal generators are: 

ve(t.x,e) " 

The vector space of infinitesimal generators associated to a translation symmetries group is denoted by T 
in the sequel. 



3.2 Constraint on infinitesimal generators of scale and translation symmetries 
— Lie symmetry determining equation 

In order to compute derivations (|18|l and thus, the associated scale symmetries group, we use the classical 
property that a derivation C associated to a symmetry of a differential system commutes with the induced 
derivation T) (in that case, the derivation C is called a symmetry of T) by extension). This fact can be state 
using Lie bracket as follows: 

[C,V]:^CoV-VoC^\V, (21) 

the parameters A is a constant; it is different from if considered symmetry does act on times and is 
otherwise. 

Infinitesimal conditions defining symmetries groups. For reader's convenience, we derive from (|21|l 
infinitesimal conditions that a derivation is an infinitesimal generators of a scale (resp. translation) symmetry 
(see § 2.5 in IL for a presentation of these infinitesimal conditions based on jet space and prolongation and |1] 
for resolution algorithms based on them). 

Lemma 1 — Infinitesimal conditions of order that an infinitesimal generator given by X]jiG(t x e) ckaJ/^/^J/ 
defines a scale symmetry of ordinary differential system E are: 

df 
y^ ayV-K^ + {at - ax,)fi = 0, for i = 1, . . . ,7i. (22) 

yG(t,x,e) ^ 

The infinitesimal conditions of order that an infinitesimal generator X^weft x B) '^v^ I^V defines a transla- 
tion symmetry of ordinary differential system E are: 

Y. «y|^=0, forz = l,...,n. (23) 

Sketch of proof. Recall that we work in vector space Der(/C/K) , thus one could consider each component of 
relation H21|) on element of canonical base (|13|l . Thus, coefficient of element d/dt is SVt — VSt = XDt and we 
deduce that the equality \ = at holds. Then, remark that the n other element's coefficient [S,Ty\Xi — XVxi 
of canonical base could be expanded as: 

(^xifi - Y "y2/"Q^ = -^fi^ for i = 1, . . . , n. (24) 

ve(t,x.e) y 

Infinitesimal conditions of order are obtained by replacing A by at in these relations. The proof of the 
second assertion is similar. D 

Example |2] (revisited) — These infinitesimal conditions associated to example 2 are given by the ma- 
tricial relation MA = with matrix M and vector A defined in figure ^ 

A:= {at a„ ap ar at^ a^a ah as a^ ) 

,. i^ fT^-fV" -- ^ ^-P i^ -^ r^\ (25) 

'lo ^ -Mo 0-^1-^0/ 

\ n n n n / 

Figure 1: Matrix defining infinitesimal conditions associated to example 2 



Exponents vector space. Considered as relations with coefficients in K{t,X,Q), the above 0th order 
infinitesimal conditions are not apparently sufficient to define completely the vector space S (resp. T) of 
exponents associated to derivation vector space S (resp. T) presented in Section l3. II and thus to determine 
searched symmetries groups. 

In fact, there is {n + £ + I) unknowns and n relations; hence there is apparently no enough relations to 
define a basis of S. 

Remark 7 — Nevertheless, this vector space is well defined because derivations V and S could be 
prolongated by the derivation 2?oo defined by formula (fTB|) (see remark IJ) and by Soo defined as follow: 

'5oo=5+^^(a,-,a,)y«|^. (26) 

jeN* yex ^ 

These derivations act on higher order derivatives of initial states variables i.e. on K{t, X, 6); The first part 
of Soc is related to relationship between classical twisted Euler derivations H18|) and scaling; the second one 
is related to the fact that 

..(^)=A(---^^,.eN,. = l,...,n. (27) 

Applying manipulations of lemmaQ]on prolongated derivations, we obtain enough relations to define desired 
vector space S by considering prolongated base elements d^ /dxi . To be more precise, the following higher 
order infinitesimal conditions hold: 

Y. "^yy—Q^ - ("-' " J'^O^'a;, = 0, (28) 

for i in (1, . . . , n) and j in N. From this infinite set of relations, one can deduce bases of S when this vector 
space is not reduced to (this is generically the case, but as shown in Section [TTTl we considered here non 
generic differential systems coming from biology, physics, etc). 

The dimension m of S gives the number of parameters of our scales symmetries group; once a basis of S is 
chosen, each of its vector is associated to a one parameter group of scale symmetries and vector's coefficients 
allows to determine exponents — the as in 117|l — of this group. 

One can choose m of the components of the ai quite arbitrarily. We shall actually want this arbitrariness 
to bear on the components corresponding to (see Section ^21). 

Remark 8 — All forthcoming computations are based and devoted to this vector space but we are going to 
see that there is no need to derive supplementary infinitesimal conditions from prolongated derivations p6|l 
and ^. 

4 Algorithm 

4.1 Infinitesimal generators computation 

Let us recall that infinitesimal conditions presented in lemma ^ show that computation of considered sym- 
metries groups is associated to computation of a kernel of a matrix of size tix (n + ^+l) with coefficients 
in the field K(t, X, <d) (this kernel defines vector space S presented at end of Section rO|l . 

Remark 9 — Fortunately, the vector space S is defined over a constant field and thus its computation 
does not require computations in K(t, X, 0). Coordinates i, X and O can be specialized to some generic 
values of a constant field in the considered matrix and so, computations could be performed numerically 
with high probability of success (i.e. computation failed when specialization are done in a Zariski closed set). 



Multiple specializations vs. higher order computations. In order to compute a basis of vector 
space S as sketched in remark |Z| one could specialize coordinates i, X, Q and their higher order derivations 
in formula H28|l i.e. compute numerical a single power series solution of linear variational system derived 
from E (see appendix I A . 1 1 for more details and ^3] for another application of that principle and more details 
on this computational strategy). This process could be done several time at different order. 

But instead, one could also perform \{n + i+ l)/n] specializations of coordinates t,X,Q in 0th order 
infinitesimal conditions presented in lemma (^ and thus, obtain a square system that allows to compute 
a base of S (see appendix I A. 2 1 for an example). Thus, we consider \{n + £ + l)/n] series computations at 
order in the sequel. 

Example [21 (continued) — After 6 specialization of matrix M defined in figure ^ one can compute 
numerically its kernel and construct the following matrix: 

/ 1 1 1 1 \ 
A' := -1 10 10 10, (29) 

\ 00-1001100/ 

such that MK = 0. Rows of this matrix define the vector space S in K^ if this vector space is given by 
coordinates (at, a„, ap, a^, cKfej, afej, afi,as, ag). By Section |31 this is sufficient to retrieve the searched 3 
parameters group (UJ of scale symmetries. In fact, the last row of matrix (|29|l is associated to derivation: 

that defines the one-parameter group of scale symmetries: p — > p/i', h -^ vh, fc2 -^ vk2- 

Computation of Gradients. Infinitesimal conditions are based on gradient computation of functions F 
defining our input differential system; we show in this section that these computations could be efficiently 
performed in our model of complexity which is described below. 

Definition 1 — Let A be a finite set of variables. A straight-line program over k[A\ is a finite sequence 
of assignments bi ^ b' Oj b" s.t. Oj is in {+, — , x , -=-} and {&', b"} is in [Jj^i{bj} U ^ U fc. Its complexity of 
evaluation is measured by its length L, which is the number of its arithm,etic operations. Hereafter, we use 
the abbreviation SLP for straight-line program. 

A SLP representing a rational expression / is a program that computes the value of / from any values of 
the ground field such that every division of the program is possible. The following constructive results taken 
form 121 allows us to determine a SLP representing the gradient of /. 

Theorem 2 — Let us consider a straight-line program computing the value of a rational expression / in 
a point of the ground field and let us denote by Lf its complexity of evaluation. One can construct a SLP of 
length 5L/ that computes the value of / and of its gradient grad(/). 

Complexity result. We gather all elements presented above and used in our computation of symmetries 
in the following proposition: 

Proposition 1 — Let E be a differential system as described in Section|2 There exists an algorithm that 
determines a rn-parameters group of scale (resp. translation) symmetries of E. The arithmetic complexity 
of this algorithm is bounded by 

o({n + t+l){L+{n + £+ l)(2n + 1+ 1))) . (31) 

Sketch of proof. Using theorem |21 the evaluation complexity of linear system derived from infinitesimal 
conditions H22(l is bounded by 0{{L + {n + (. + \)n)n). This system should be evaluated \{n + (. + l)/?^] 



times on some generic specializations elements in a constant field in order to obtain a square system. Using 
Gaussian elimination, a base of the resulting system kernel could be computed with 0[{n+i+l)^) arithmetic 
operation. Thus, we obtain our complexity l|31|) . D 

Remark 10 — The algorithm presented in this section is not probabilistic even if it is based on special- 
izations that could lead to computation of some spurious symmetries. In fact, it could occurs — with a very 
small probability — that coordinates t, X, & are specialized on several points in the orbit of a symmetries 
group of S; in this case, some spurious symmetries are obtain. 

Fortunately, a simple evaluation of computed symmetries on our original system allows to show if it is a 
computational artifact or not. The complexity associated to these tests is bounded by 0{{n + £ + 1)L). 

Remark 11 — There is an infinite way to choose a basis of S. But, one can use Lenstra, Lenstra and 
Lovasz' basis reduction algorithm in order to obtain a reduced basis in the sense that exponent as are smaller 
then what could be obtained using classical Gram-Schmidt orthogonalization. 

4.2 Computation of some rational invariants and original system rewritting 

Generally speaking, when a differential system admits a group action as symmetry, we can rewrite it in 
terms of the invariants of the group action. Fels and Olver revised the moving frame the construction for 
(differential) invariants ^Jj. The constructed invariants allow for a trivial rewriting. We do not recall the 
general theory but rather show how it works on the group actions of interest in this paper. 

While computation of non-specific rational invariants i.e. a total description of invariant field for any al- 
gebraic group action could be done using reduced Grobner basis computation (see j7]), we do not need such a 
general tool because we restrict ourself to scale and translation transformations. In fact for scale/translation 
symmetries, we could restrict ourself to compute £ — m time independent rational invariants in the multi- 
plicative group in K(8) generated by the set: 

M := {6"^ I (6*,/?) c e X Q, V5 e S, S0 ^ O}. (32) 

A set of generators of this multiplicative group is denoted by tti, . . . , Tr^^m- Furthermore, we restrict ourself 
to looking for following kind of time dependent rational invariants: 

TTt = ptt, TT^. = p^Xi, I <i <n, (33) 

where the pi are in M. This arbitrary choice simplify invariants' computation and is well suited with our 
purpose; we are trying to use computed group of symmetries to reduce the number of parameters and not 
to reduce the number of equations. 

Group action. Given o'i\i....^\^^) a ?n-parameters symmetries group, let us consider its action 'ijj on 
fieldK(t,X,e): 

F™ xK(t,x,e) ^ K(i,x,e), 

(Ai,...,A,„) X y -^ o-(Ai,...,A„)(y)- 

with F a field where parameters could be specified. For scale symmetries, we could wrote ^ as the substitution 
given by: 

Vye(i,X,e), V^(y)=2/niLiAf-, (35) 

where exponents are given by the basis {{at,i, ax^.i, ■ ■ ■ , clxiA, cloiAi ■ ■ ■ , aei.i), i = 1, . . . , m\ of S. 



10 



Rational invariant computation. By classical canonical homomorphism, the multiplicative set F™ x M 
could be considered as a Z module of dimension 2€. Thus, the subset of relations 1)35(1 involving only 
parameters could be represented by the following £ x 2£ matrix: 





Ai 


Am 


Oi 




( ae,,i 




aei,m 


1 





092,1 




0.02, m 





1 


oe,_i,i 




0'9i_i,m 






\ «e,,i 




0,01, m 








\ 



1 
1/ 

One can consider the matrix obtained after a permutation of 1(3611 lines : 

Ai ... Am 01 ... 9e 



(36) 



/ 



as. 



aa 



V 



1 
1 







\ 



(37) 



such in order to ensure that the determinant of the submatrix (a^ -Yi^i 



is not 0. 
A Gaussian elimination performed on this matrix and terminated at to + 1 column position leads to the 
matrix : 

7i,ei •■• 7i,e, \ 



/ 1 









Pm+l,e: 






(38) 



V ■ ■ • f^iA ■ ■ ■ (^iA I 
This computation is sufficient to determine the following generators of the multiplicative set ((32(1 of rational 
invariants: 



'^(Ai 




>o-(Ai,...,A„) 




(39) 



Thus, this elimination construct elements tt^ of /C that are invariant under the action -0; remark that con- 
sidering the whole action of -0 i.e. a (n + ^ + 1) x 2(n + ^ + 1) matrix, the same process allows to construct 
time dependent invariants of type ((33(1 but we are not going to do so. 

Example 121 (continued) — In Section ITTl we determine by numerical computation, a basis of vector 
space S presented in matrix 1(29(1 : taking the 6 last columns of this matrix, one obtains the 3 first columns 
of the following matrix: 

/010100000\ 

10 10 

110 10 

10 10 ■ (40) 

10 10 

yioooooooiy 
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After a Gaussian elimination on the first 3 columns of above matrix, we obtain: 



(41) 



/OOOl 0-1 100\ 

10 1 

10 1-10 

10 10 

0-1 110 
\0000-l 001/ 

The first line of above matrix defines the rational invariants f) = rh/k2 and the last line defines e = e/ki 
etc. 



Algebraic counterpart of Probenius theorem in considered cases. First, let us remark that the 
action ip is surjectivc. In fact, one can consider the group's parameters Ai = . . . = A„j = exp(O) = 1 that 
correspond to identity map. This map is associated to derivation in S by formula (|19|l . 

The arbitrary choice of parameters made at beginning of Section 14.21 is motivated by the possibility to 
rewrite our original dynamic Hll|l in order to reduce the number of parameters; we want to determine an 
expression of the dynamic V on an intermediate invariant field T = Kf i.e. 



T = ^{■KnTrxi, ■ ■ ■ -.T^x^.T^X-, ■ ■ ■ ,T^e-m) 



JC, 



which transcendence dimension w.r.t. K is smaller then the original one but where the number of variables 
depending on times does not change. 

In the following proposition, we gives an elementary algebraic proof of a general results for reader's 
convenience (see 7 for definition, computation and use of replacement invariants). 



Proposition 2 



There exists Ai, . . . , Am in M such that the following equality holds : 



i^-Hi^"^'' ^".)) = (ai,...,a„) 



x/C. 



(42) 



Sketch of proof. This proposition is a reformulation of matrix (|38|l 's structure inherited from matrix (|36(l . 

YiZi' 'm keep a track of the performed gaussian elimination i.e. the following 



In fact, the submatrix (7. 
equality holds: 



')'iii 



Ti.Sj 



/I 




as 



If, for i — 1, . . . ,m 
relations hold: 



ax 







: 

1 



o\ 



(43) 



V ... 1 ... / 

we define A,; := Y[j=i ^j ' ^ ^^^ notice that the above equality shows that the following 



cr. 



'(Ai,...,A„) 

Therefore, there exists a subset 
ions 



)-^n H' 



'^h,e, 



^n^ 



Z^ h. = l '^6i,h'^h,i 



(44) 



h=i \i=i 



i=i 



the relations a,^ ^ AOi) = 1 hold for i = 1 



6rn) of parameter's set and a subset (Ai, . . . , Am) of K such that 
,m. The same type of result is valid for the subma- 



trix(A,,)^-^;;t..,,ofeHl 

(ax 



ae„+i,mV^i.«i 



V "e 



\7„ 






■ ■ l-'3„+l,«„ 



(45) 



^'f^ij, J 
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This matricial relation prove that, for i = rn + 1, . . . , £ the foUowing equahties hold: 



'^(Ax 







To conclude, remark that the same properties hold for time dependent variables. D 

Thus, after Gaussian elimination performed on H36() . we obtain a description of /C^ and an applica- 
tion (T,c c , that maps K- to IC^ . In fact, This action allows to determine rational invariants and rewrite 
original system in an set of invariant coordinates with a reduce number of parameters. 

Example 121 (continued) — The second, third and fourth lines of matrix 141|l give the specialization ||SJ) 
of parameters that allows to determine time dependent invariants and system rewriting. We summarized 
computation done in this section by the following proposition: 

Proposition 3 — Computation of rational invariants defined in (|32|) and by (|33ll could be performed by 
Gaussian elimination and thus its complexity is bounded by 0[£^). 

Remark 12 — The computation of rational invariants presented in this section suppose that the considered 
symmetries groups is composed of dilatation. The same consideration holds for translation. In fact, the group 
action is in this case 

n 

Vye(t,X,e), i;{y)^y + J2ay,^X^, (47) 

and is thus linear; that allows exactly the same type of computation. 

Remark 13 — If S and T are two symmetries of V, the Lie bracket [S, T] is also a symmetry of V by 
Jacobi identity. Thus, symmetries form a Lie algebra and if Si is in S and % is in T, we have the classical 
facts that [5i, 5j] = and [7^, Tj] and [5^,7^] are in T i.e. induced derivation is of the same type then (|2U|) . 
Typically, scaling are symmetries of translation's invariant while the opposite is not true. 

Thus, S and T generate a solvable Lie algebra and above commutation relations show that we have to use 
translation symmetries first in our algorithm to reduce parameter's number and then use scale symmetries 
(see § 2.5 in IT). 

Remark 14 — As mentioned in remark [TTl there is some freedom in choosing the components of a basis 
of S (resp. T), the set of exponents that define a scaling (resp. translation) symmetry of a differential system. 
We shall always try to have the freedom on the components of a that correspond to 8. This is achieved 
by placing correctly the unknown when solving the system by Gaussian elimination; this fact allows user 
to choose parameters to eliminate. When this is not achievable, the general Lie method's could be use to 
solve — partially — the system by quadrature. 



5 Conclusion and future work 

In this paper, we consider the computation of scale and/or transformation group that are symmetries of 
an ordinary differential system and the determination of some of their invariants. We use these groups and 
associated invariants in order to rewrite the ordinary differential system in an set of invariant coordinates 
with less parameters. The complexity of this process is polynomial in input's size. 

Remark 15 — For the sake of simplicity, we do not include control time dependent variables U in our input 
system (|11() : if such variables occur, the ground field is in practice the differential fraction field K(t/). Com- 
putations are performed after specialization of variables U on power series with random integer coefficients 
and which are truncated at order n + £ + 1. 
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Same type of result could be proved for more general symmetries that allow to reduce further the number 
of significant parameters as shown below. 

Example 3 — Let us consider a FitzHugh Nagumo model (see § 7 in ^U]) defined as follow: 

a = b = c = d = 0, dx/dt = c{x-x^/3-y + d), dy/dt = {x + a - by)/c. (48) 

This system does not have scale or translation symmetries that are considered in this paper but one can 
determine that the derivation: 

— h— — 
dy da dd 

is an infinitesimal generators of the one-parameter symmetries group (which is not of type H2U|) '): 

y -^ y + X, a ^ a + bX, d -^ d + X. (50) 

Up to our knowledge, there is likely no polynomial time algorithm that compute infinitesimal generators (|49(l . 
In fact, this type of symmetries can be found by supposing that all coefficients of seeked infinitesimal 
generators of symetries are rational function of parameters; in that case the solution of determining system 
of PDE H21|l is reduced to the computation of a polynomial matrix kernel. But as done in this paper, one 
can use the following invariant coordinates t) ^ y — d and a ~ a + bd to rewrite system ()48f) as follow: 

dx/dt = c{x-x^/3~t)), dt)/dt= {x + a-bt))/c (51) 



Acknowledgments. The second author is grateful to F. Lemaire and M. Safey El Din for many useful 
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A Two computation methods 

Considering a classical example taken from biology, we are going to explicit some computation evoked in 
Section |4. II First, we determine a symmetries group using power series approach and then, we retrieve the 
same result using just specialization of 0th order infinitesimal conditions. 

A.l Series based computations 

Let us consider the linear variational system: 

S = F{t,^,Q), 
d 5S dF , „ ^^ dE 

dtdx = ax^^'^'^fe' (52) 

d dE OF , „ ^,dE dF , „ ^, 

with the initial conditions dE/dX = Id„xri and dE/dQ = when t = 0. Power series solutions of this system 
are (I15|l and: 

dE _^ dV^X P dE _y-^ dV^X P 

dX ~ -^ dX 71' de~ ^ dQ j! ' ^ ' 

Coefficients of these series are used in generalized infinitesimal conditions H28|l that allows symmetries com- 
putations. 

Hence, one can construct a system of ordinary differential equations that allows to compute directly a 
specialization of T>^ Xi and dV^Xi/dy with y in {X, Q). In fact, on can compute power series solutions of this 
system and, the wanted quantities are coefficients of these power series. 

Example 4 — We perform above computations on the following Michaelis Menten's equation (see § 6.3 
in dm): 

i^T^ = f{ki,k2,0- (54) 
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As df /dt = 0, the associated linear variational system is: 



d 9« 
dt dx 

d 94 


kik2 di 1 C 


dt 9fci 

d 94 
dt 9fc2 


(fc2+?)^afcl ' (fc2+«)' 

(fe2+S)^ 9fc2 (fc2+e)^ 



(55) 



Using specializations defined by ki = 7 , k2 — 2 and ^(0) = 3, at order 5 power series solutions ^(i, 3, 7, 2) of 
above equations are: 



147+2 



1372 43 



<; "^ 4" 5 '' ^ 125 '' 3125'' 



.^L — 3j- , 42 .2 
9fei ~ 5 ^ 125 



588 J.3 



3125 



t^ 



15625 



2401 ^4 , (5/^5^ 
31250'' ^ ^\^ h 

t^ + Oit^), 



9; 

9fc2 

ax 



25 

14 



t 



1250 ^ 
196^2 
625 



1029 ^3 



69629 +4 



3125 
I _68l 
^ 9375 



312500 



es^t^ + ^^^t^ + Oit^). 



(56) 



25 625 ' 9375 ' 234375 

Using these values, we construct linear system associated to commutation condition l|21|) : 



/ 



\- 



-1 







63 

25 





21 



\ 



294 
125 

8232 4116 
3125 3125 



2646 588 
625 125 



_42 
25 

294 
'625 



16464 12348 
3125 3125 



fa,\ ( 
= \ 






1 

21 

294 
125 



\ 



(57) 



\Oik2/ \ 8232 / 



and by a Gaussian elimination, we determine the basis of S defined by the relations at — — A,Q!fej ~ X + fi 
and ak2 = fJ-iCtx — fJ-- Thus, the dimension of vector space S is 2 and the following derivations: 



„ 9,9 „ d , d ^ d 



(58) 



form one of its bases (for the sake of simplicity, higher order derivatives are omitted in these derivations). 
We deduce from p9f) that the following 2 parameters groups: 



t 

X 



\t, 
HX, 



fci 

fc2 



^ki/X, 



(59) 



acts on (|54|l and that t = kit/k2 and y = x/k2 are a convenient set of new coordinates in which system (|5 
could be rewritten as dy/dt = j:/(l + y). 



Using specialization: 



A. 2 Multiple specialization 

Example 0] (continued) — For this example, 0th order infinitesimal condition (|22|l leads to consider the 
following vector: 

( k\x fci x'^ fci X xki ^2 \ , , 

\k2+x'~ [k2+xf' k2+x'~ {k2+xf) ' 

X = -2, ki = 10, fc2 = -2, x = -4, ki = -7, k2 = 1, 
a; = 2, /ci = 8, ^2 = -1, x = A,ki = -2, fc2 = 1 

one can obtain the following matrix associated to infinitesimal condition H22|l : 

/ 5 -5/2 5 -5/2 \ 

-28/3 112/9 -28/3 -28/9 

16 -32 16 16 

V -8/5 32/25 -8/5 8/25 / 



(61) 
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whose kernel is defined by vectors (1,1,0,1) and (—1,0,1,0) that generates the vector space S already 
determined above (applying LLL reduction on these vectors, we retrieve exactly the previous basis). 

Computation of numerical power series solutions of (|54|l and H55|l performed above could be done using 
specialization H61() . Instead of considered series solutions associated to a single specialization, one can 
consider two or more such series. The series' order needed by our computation decreases with the number 
of used specialization. 

So, there is no need to compute power series as described above even if theoretical structures are clearly 
defined using this approach (see remark [TJ. 
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